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Observations of neutrino oscillations show that neutrinos have mass. However, the best constraints 
on this mass currently come from cosmology, via measurements of the cosmic microwave background 
and large scale structure. In this paper, we explore the prospects for using low-frequency radio 
observations of the redshifted 21 cm signal from the epoch of reionization to further constrain 
neutrino masses. We use the Fisher matrix formalism to compare future galaxy surveys and 21 cm 
experiments. We show that by pushing to smaller scales and probing a considerably larger volume 
21 cm experiments can provide stronger constraints on neutrino masses than even very large galaxy 
surveys. Finally, we consider the possibility of going beyond measurements of the total neutrino mass 
to constraining the mass hierarchies. For a futuristic, 21 cm experiment we show that individual 
neutrino masses could be measured separately from the total neutrino mass. 
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I. INTRODUCTION 

The potential for using cosmological perturbations to 
constrain particle physics in general and neutrinos in 
particular was pointed out more than twenty years ago 
[ij. Recent cosmological data have indeed allowed much 
progress in constraining dark matter properties (e.g. 
0, [| ) and neutrino physics ( 0, i, i, 0, H i,M [HI, [ll ) ■ 

Perhaps the most interesting neutrino property that 
can be investigated with cosmology is its mass. While so- 
lar and atmospheric neutrino oscillation experiments set 
a milestone for particle physics models by showing that 
neutrinos do have a mass, by measuring mass differences 
(Amii ~ 8 X lO-^eV^ and Amii ^2x IQ-^eV^; [H) 
they don't carry any information on the absolute neu- 
trino mass scale. Interesting limits on the neutrino mass 
scale may come from tritium beta decay in the near fu- 
ture (e.g. KATRIN), which will probe electron neutrino 
masses above 0.2 eV. 

Recent cosmological data, however, also provide the 
opportunity to probe neutrino masses through the pecu- 
liar behavior of massive neutrinos, which are relativistic 
for part of the Universe's history and possibly not at 
later times. This behavior has an impact on both the 
background evolution and the growth of structures on 
intermediate and small scales. The constraints are typi- 
cally derived by combining information on the amplitude 
of cosmological perturbations at the redshift of recombi- 
nation (probed by the Cosmic Microwave Background) 
and in the very local Universe (probed by galaxy surveys 
such as the 2dF or SDSS). Cosmology currently sets the 
tightest upper hmit on the neutrino mass (e.g. [1, 
to be about 0.62 eV. It should be noted here that this 
same mass limit was recently derived with the WMAP- 
5 year data combined with limits from supernovae (SN) 
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and baryon acoustic oscillations (BAG) With fu- 

ture surveys, alternative techniques will be available to 
determine neutrino mass, including the use of CMB-LSS 
cross correlation fl^ and lensing studies [T7| . 

While most methods mentioned above make use of the 
fluctuation amplitude either at high redshifts (z ~ 1000) 
or at low ones {z ~ 1), a new cosmological observable 
has emerged in the past few years as a potential probe 
of fluctuations at intermediate redshifts (6 < z < 20): 
the study of fluctuations in the 21 cm line. As the power 
spectrum of perturbation at these redshifts is affected by 
the presence of massive neutrinos, the study of the 21 
cm line may provide complementary information on neu- 
trino masses [l^. Moreover, as fluctuations are smaller 
at higher redshifts, the effects of non-linearity at small 
scales are weaker than at the present time, allowing for a 
more straightforward interpretation of the data and ulti- 
mately to the use of a broader scale range. 

As the Universe is still matter dominated at such high 
redshifts, cosmological probes of neutrino masses with 
the study of the 21 cm line also offer the advantage of 
avoiding the typical degeneracy between neutrino mass 
and dark energy which is seen in probes relyin g on power 
spectrum measurements at redshifts z < 1 [19|. The 
study of 21 cm fluctuation would, however, still exploit 
the impact of neutrinos on the matter power spectrum 
by probing perturbation amplitudes at redshifts substan- 
tially different from the ones probed by CMB studies. 

The cosmological limits cited above were derived under 
the hypothesis that neutrinos all share the same mass. As 
neutrino oscillations only determine two of the three mass 
square differences, we do not know the precise hierarchy 
for the masses (toi < m2 < ma or < mi < m2 repre- 
senting normal and inverted hierarchy respectively) . The 
relation between the mass of the lightest neutrino and the 
sum of all neutrino masses clearly depends on the par- 
ticular hierarchy, however this relation is affected by the 
hierarchy only if the lightest neutrino has a mass below 
~ 0.3eV. Also, for the time being, the effect of different 
neutrino masses has been largely irrelevant for cosmo- 
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logical considerations as the upper limits on the mass 
imposed by cosmology are well above the mass splitting. 
As cosmological data become more and more abundant 
and different probes contribute in determining the neu- 
trino mass, it is possible that cosmology will be sensitive 
to the masses, which are comparable to the measured 
mass splitting. In such a case, cosmology would, not 
only, set the tightest measurement of the neutrino mass 
scale, but could also enable one to distinguish between 
the normal and inverted hierarchy (hereafter NH and IH 
respectively) . 

In this paper, we consider the potential of constrain- 
ing neutrino masses with future surveys mapping the 21 
cm line. We will investigate how these measurements 
complement and improve the results based on the future 
CMB probes (Planck) and large scale structure surveys. 
In particular, as previous studies on this topic suggest 
that the precision in the determination of the neutrino 
mass should be lower than the mass splitting [20[, we 
study the potential of future 21 cm surveys to determine 
the neutrino hierarchy. 

This paper is organized as follows: in sec. |TT] we dis- 
cuss the generalities of the 21 cm probe; in sec. IIIII we 
discuss how the neutrino mass specifications affect the 
power spectrum at the relevant scales and redshifts, in 
sec. II VI we present the adopted Fisher matrix approach, 
while sec. IVlandlVIIare dedicated to results and conclu- 
sions respectively. 



II. 21 CM AS A COSMOLOGICAL PROBE 

We begin by briefly summarizing the physics of the 21 
cm signal and refer the interested reader to Furlanetto, 
Oh, and Briggs [2l| for further details. The 21 cm line of 
the hydrogen atom results from hyperfine splitting of the 
IS ground state due to the interaction of the magnetic 
moments of the proton and the electron. The HI spin 
temperature Ts is defined through the ratio between the 
number densities of hydrogen atoms in the IS* triplet and 
15 singlet levels, ni/no = {gi/ go) e'xjp{—T-i,/Ts), where 
{91/90) — 3 is the ratio of the spin degeneracy factors 
of the two levels, and = /ic/fcAaicm = 0.068 K. The 
optical depth of this transition is small at all relevant 
redshifts, yielding brightness temperature fluctuations 



n = 27XHI ( 1 + g'^b 
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Here xhi is the neutral fraction of hydrogen, 6b is the 
fractional overdensity in baryons, and the factor of 4/3 
arises from including the effect of peculiar velocities and 
averaging the signal over angles. Provided that there is a 
sufficient background of Lya photons, the spin tempera- 
ture Ts is coupled to the gas kinetic temperature Tk via 
the Wouthysen-Field effect [22,[2i|. 



In general, Tf, is sensitive to details of the IGM and 
is complex to model. However, with a few reasonable 
assumptions this dependence drops out. At redshifts 
long after star formation begins, it is likely that the X- 
ray background produced by early stellar remnants has 
heated the inter-galactic medium (IGM) so that the gas 
kinetic temperature Tk is much greater than the CMB 
temperature Tj [U, [2^, [2^. The same star formation 
produces a large background of Lya photons sufficient to 
couple Ts to Tk- In this scenario, appropriate for z < 10, 
we are justified in taking Ts ~ Tk 3> Tj in which limit 
Tfc shows no dependence on Ts- We need then only worry 
about fluctuations in the density and neutral fraction. If 
the IGM is fully neutral, 21 cm fluctuations arise only 
from the variation in density from point to point and P21 
provides a tracer of the matter power spectrum. In this 
limit, we can write the power spectrum of 21 cm bright- 
ness fluctuations as 



P2i{k) = T^Pss{k){l 



2\2 



(2) 



where Pss is the matter power spectrum and /x is the an- 
gle between the line of sight and the wavenumber k. In 
this situation, 21 cm observations may be thought of as 
being analogous to galaxy surveys with Tb replacing the 
unknown galaxy bias. While in principle Tb can be calcu- 
lated given a knowledge of cosmological parameters, its 
dependence on the largely unknown ionization and ther- 
mal history means that we should not treat it as a known 
quantity, but instead as a free parameter to be measured. 
Thus cosmological information will come primarily from 
the shape of P21 and not its amplitude. 

Throughout this paper, we will adopt the above op- 
timistic assumption that 21 cm observations at 2; < 10 
measure the underlying matter power spectrum modu- 
lated by the 21 cm brightness temperature, i.e. P2i{k) cx 
T^P{k). In practice, this is unlikely to be the case since 
a combination of fluctuations in the Lya flux [13, HI], 
inhomogeneous X-ray heating [2^ . and inhomogeneous 
reionization ^SflJ will modify the 21 cm power spectrum 
signiflcantly. Although these fluctuations are unlikely to 
all contribute at the same time [26l|, finding a redshift 
range where only the density fluctuations are significant 
would be extremely fortuitous. In principle, the matter 
power spectrum may be measured in the presence of these 
other contributions by precise measurements of the angu- 
lar dependence of the 21 cm power spectrum [3l[. This 
procedure, however, reduces the information extracted, 
therefore the constraints on cosmological parameters that 
may be achieved in this pessimistic scenario are seriously 
weakened [S^l • More likely detailed modeling of reioniza- 
tion will allow cosmological parameters to be extracted 
in the presence of neutral fraction fluctuations with only 
moderate degradation ~ 10 — 50% in the constraints [33 |. 
Finally, since the majority of these extra sources of fluctu- 
ation are most signiflcant on large scales, on small scales 
k > 1 Mpc"^ it might be hoped P2i{k) will trace P{k) 
closely (although see Lidz et al. 33]). 

We note that, although the details of reionization are 
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still poorly constrained, increasingly accurate observa- 
tions of the CMB constrain the optical depth to the sur- 
face of last scattering. This provides an integral con- 
straint on the ionization history that can be used to 
limit when reionization ends. The most recent WMAP 
results [3| give r = 0.09 ± 0.017. If reionization is as- 
sumed to end instantaneously this translates into Zroion = 
10.8 ± 1.4, while more extended reionization models lead 
to reionization completing at lower redshifts (e.g. us- 
ing the model of Furlanetto [1^ gives Zreion = 9 ± 1.4). 
A complementary constraint on the end of reionization 
comes from observations of the Gunn-Peterson trough in 
the Lya forest in absorption spectra of quasars at red- 
shift z > 6.5 [s^, [sEl- This is widely interpreted as in- 
dicating a neutral fraction of xhi ^ 10~^ in the IGM, 
possibly signalling the tail end of reionization. Guided 
by these observations, we will model reionization as oc- 
curring sharply at z = 7.5. Although simplistic, this sets 
a reasonable lower limit on 21 cm observations during 
the epoch of reionization. Although some neutral hy- 
drogen remains in dense regions after reionization and 
may be targeted by future 21 cm observations [sHi, the 
post-reionization signal strength decreases significantly. 
For this reason, we will restrict our 21 cm experiments 
observations to z > 7.5. The upper limit on the acces- 
sible redshift range comes from the experimental design. 
For currently proposed experiments, this is likely to be 
around 100 MHz (z = 13). Given the rapidly increas- 
ing foregrounds and uncertainty in the spin temperature 
at higher redshifts, we will limit ourselves to redshifts 
z < 10. In principle, increasing the redshift range will 
increase parameter constraints, but the increasing ther- 
mal noise reduces the utility of additional high redshift 
bins. 



III. EFFECT ON NEUTRINO MASSES ON THE 
POWER SPECTRUM 

Although the physics of neutrinos is potentially very 
rich and models exist for which the number of neutrino 
species is different from three, for the purpose of this 
paper, we will only consider the standard case with three 
light neutrino species, and will discuss the implications 
of having different assortments of neutrino masses. 

If neutrinos were in thermal contact with the CMB 
at early times, their number density is fixed. Therefore, 
their contribution to the total present energy density is 
strictly related to their masses: Vli, — mi/93. 14ft,^eV 
and they would contribute a fraction = ^l^/^lm to 
the total mass density of the Universe. From neutrino 
oscillations, we know that at least two neutrinos are 
non-relativistic today, as the mass differences implied are 
about 0.047 and 0.009 eV. The current Hmit on the neu- 
trino mass scale is about 0.6 eV. If neutrinos have masses 
above 1.7 x 10~^eV (that is the neutrino temperature to- 
day) and below the current mass scale limit, they are 
still relativistic when cosmologically relevant scales enter 



the horizon. Therefore, because of their free-streaming 
behavior, they then tend to erase fluctuations on such 
scales. For a given neutrino mass, the impact of the free- 
streaming behavior on the actual matter power spectrum 
depends on the scale considered: scales larger than the 
horizon size when neutrinos become non-relativistic are 
unaffected; while smaller scales suffer a damping. As- 
suming that neutrinos contribute to the expansion of the 
Universe, but not to the clustering, the growth of cold 
dark matter density perturbations as a function of time 
is: ScDM oc ai-(3/5)/„_ pgj. gjj^all values of the neutrino 
mass, this effect is tiny; However, it may last for a very 
extended period of time, resulting in a sizable decrement 
of the matter power spectrum. Since the CMB probes 
a perturbation's amplitude at decoupling (z ~ 1000) 
and the 21 cm measurements probe much lower redshifts 
(z ~ 6 — 10) the effect can become significant. Free- 
streaming is not the only physical effect caused by mas- 
sive neutrinos impacting the matter power spectrum. As 
neutrinos are relativistic at early times and become non- 
relativistic at late times, they contribute to the radiation 
energy density initially and to the matter energy density 
at later times. Therefore, according to the actual values 
of their masses, they may slightly change the redshift of 
matter-radiation equivalence; which in turn has an im- 
pact on when dark matter perturbations start growing at 
the rate of a matter dominated Universe. Fig. [1] shows 
the effect of neutrino mass on the power spectrum at 
different scales and redshifts, by looking at the deriva- 
tive of the power spectrum with respect to M^, assuming 
the three species have total mass = 0.3 eV. Smaller 
scales are the most affected by different neutrino masses. 
Such scales tend to be non-linear at the present time, 
but still linear at z = 8. The 21 cm probe can therefore 
take advantage from the information provided by all such 
small scales. 

Apart from the general effects implied by neutrinos 
having a mass, should the overall mass scale be below 
~ 0.3 eV, then the behavior of each single neutrino 
species may impact the power spectrum. This hypoth- 
esis has been investigated in detail in [37| . which shows 
maximum differences of the order of 0.5% in the present 
power spectrum range 0.001 < fc < 0.1 /iMpc"^ for differ- 
ent mass configurations and a total fiducial mass of 0.12 
eV. Such small differences cannot be probed with current 
galaxy surveys like SDSS, however they may be with fu- 
ture cosmic shear surveys. We argue here that future 21 
cm surveys may well also allow us to distinguish between 
different mass configurations, in particular, it may en- 
able us to distinguish between the normal and inverted 
hierarchy. Fig. [2] represents the difference in the power 
spectra between the normal and inverted hierarchy rela- 
tive to the case of three neutrinos with degenerate masses 
at redshifts z = and z = 8, of interest for the 21 cm 
probe. Although differences in the power spectrum are 
tiny (0.1 — 0.2%), at z = 8 they extend in the precise k 
range probed by the 21 cm surveys: 0.1 < k < 2 hMpc~^. 
In the following, we will discuss which kind of surveys 
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FIG. 1: Effect of a different neutrino mass on tlie power spec- 
trum at redshifts 2: = 0.3 (soiid curve), 4 (daslied curve), and 
8 (dotted curve). Tlie curves represent tlie derivatives of tlie 
matter power spectrum witfi respect to neutrino mass for a 
fiducial mass value of Mi, =0.3 eV. The vertical dashed lines 
represent the non- linear scale at redshifts z = 0.3, 4, and 8 
(from left to right). Note that the derivative becomes more 
negative at lower redshifts, as massive neutrinos have more 
time to affect the power spectrum. 



FIG. 2: Difference between the power spectra between the 
normal and inverted hierarchy at different redshifts. We plot 
the ratio of the power spectrum calculated with total neutrino 
mass Ml, =0.12 distributed in the normal hierarchy (dotted 
curve) and inverted hierarchy (solid curve) relative to the case 
where the three neutrino masses are degenerate (mi = 7712 = 
ms). These curves are calculated at z = (thin curves) and 
z = 8 (thick curves). The difference continues to smaller 
scales, but saturates at ~ 0.2% at = l/iMpc~^. 



may reach the sensitivity to neutrino masses which allow 
the determination of the neutrino mass hierarchy. 

Another advantage of studying the matter perturba- 
tion evolution through the 21 cm line is that at the rel- 
evant times the Universe is still matter dominated. At 
lower redshifts z <2, dark energy dominates the expan- 
sion and influence the perturbations' growth, slowing it 
down with respect to the matter dominated regime. The 
actual growth rate depend on the specific value of w{z)^ 
the equation of state parameter for dark energy, which is 
one of the least currently constrained parameters. As a 
result, for a given perturbation' s amplitude at the time 
of recombination, an observed value at the present time 
showing a slower evolution than that implied by a matter 
dominated model may be interpreted as either a (small 
but long-lasting) neutrino effect or a (big, late time) dark 
energy one. There is currently no way of decoupling the 
damping effect of neutrino free-streaming on the power 
spectrum from the dark energy effect, unless other probes 
not involving growth of perturbations can characterize 
the evolution of dark energy very well. 

We note here that current limits on neutrino masses 
at the level on ~ 0.6 eV are not only reached using 
galaxy surveys and the Lyman-a forest, but also using 
WMAP-5 yrs in combination with BAO and SN mea- 
surements, therefore completely avoiding the use of the 



matter power spectrum [15|]. In this case, the degener- 
acy between the neutrino mass determination and the 
quintessence equation of state no longer exists. Addi- 
tionally, 21 cm studies will be subject to a completely 
different set of systematic uncertainties to galaxy surveys 
and the Lyman-a forest. Exploiting multiple techniques 
to probe the matter power spectrum on small scales will 
likely be necessary to ensure that the various systematics 
of different techniques are under control. 

However, it is unlikely that the CMB alone will ever 
greatly improve on this constraint; leaving the fine mea- 
surements of neutrino masses to some kind of probe in- 
volving also the matter power spectrum. The 21 cm 
probe capitalizes on the effect of neutrino masses on 
small scale perturbations while avoiding the low-redshift 
complications implied by the dominance of the currently 
poorly-characterized dark energy and also the greater 
non-linearities existing on small scales. 

Let us finally comment on the role of the CMB in de- 
termining the neutrino mass. Given that the CMB de- 
couples when the temperature of the Universe was 0.26 
eV, it cannot be sensitive to very small neutrino masses, 
with the possible exception of when the lensing effect by 
lower redshift structures is taken in consideration. How- 
ever, the upcoming Planck experiment has an incredible 
constraining power on generally all cosmological param- 
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eters, which is hardy attained by any other single cosmo- 
logical probe. In turn, Planck will enable other probes 
(like galaxy surveys and the 21 cm) to attain a higher pre- 
cision in the neutrino mass determination by providing 
accurate measurements for other cosmological parame- 
ters like flrn and Hq with which the neutrino mass may 
be partially degenerate. We will therefore use Planck as 
a baseline and will consider its projected performances 
in our calculations. 



IV. FISHER CALCULATION 



TABLE I: Specification for CMB experiments. For each chan- 
nel we list the frequency in GHz, the beam size Sbcam is the 
FWHM in arcmin. Sensitivities ar and ap are in fiK per 
FWHM beam, w = (^boamO")"^. Taken from Lesgourgues 
et al. In each case, we assume /sky = 0.65 allowing for 

subtraction of the galactic plane. 



Experiment 


Frequency 


^bcam 


(JT 


ap 


Planck 


100 


9.5 


6.8 


10.9 




143 


7.1 


6.0 


11.4 




217 


5.0 


13.1 


26.7 



To make predictions for parameter constraints from fu- 
ture surveys we make use of the Fisher matrix formalism. 
This has been described in detail elsewhere [S^ and we 
provide only a brief summary here of the details relevant 
for our analysis. 

We take as our basic set of parameters the three pa- 
rameters describing the energy density in matter Vlj^h?, 
baryons 0,bh? , and dark energy SIa, two parameters 
describing the tilt and amplitude of the primordial 
inflationary power spectrum As^ and the optical depth 
T to the surface of last scattering. In addition to these 
basic parameters, we add the total mass in neutrinos 
Ml, = mi + 7712 + W3 and, since it modifies power on 
small scales the helium fraction Yhc- We assume a flat 
Universe throughout. This gives a parameter vector 
e = {Qrnh'^,^bh^,^A,ns,Al,T,Yuc,M„). In addition, 
we require an effective bias - either a genuine bias 
for galaxy surveys or the brightness temperature for 
21 cm experiments. We will separately explore the 
degeneracies introduced by allowing the equation of 
state w of the dark energy to vary. Throughout, we 
use CAMB [s^ to calculate CMB anisotropics and 
matter power spectra. Our fiducial set of parame- 
ters is e = {n,nh'^,nbh'^,nA,ns,A'^,,T,YHc,M„) = 
(0.147, 0.023, 0.7, 0.95, 26, 6, 0.1, 0.24, 0.3 eV), which 
yields as — 0.77, consistent with the latest WMAP data 

To evaluate Fij, we need to specify a model, which 
determines the dependence of the likelihood function on 
0, and a point in parameter space where we wish to 
determine parameter uncertainties. In the case that the 
model parameters are Gaussian distributed, the Fisher 
matrix takes the form 



-.. = i..(c-.c..c-.c„, + |ic-.|i, 



(3) 



where C is the covariance matrix for the data, and n is 
the data's mean. This will be a good approximation in 
the case of CMB observations, galaxy surveys (see [i^l 
for a discussion on full likelihood constraints), and 21 cm 
observations. Note that, for our purposes, we will need 
to combine information from these three sources. When 
used together these data sets break many degeneracies 
that are present when they are used alone. Let us con- 
sider the Fisher matrix from each in turn. 



A CMB experiment may be characterized by a beam 
size ^beam and sensitivities to temperature <tt and polar- 
ization ap. Given these quantities, the Fisher matrix is 
given by [ill, El S, SI 



I X.Y " " 



(4) 



where Cf^ is the power in the ^th multipole for X = 
T, E, B, and C, the temperature, E-mode polarization, 
B-mode polarization, and TE cross-correlation respec- 
tively. The elements of the covariance matrix Cov£ be- 
tween the various power spectra are given in [43l . |4^ . 
Since in the near future the best CMB observations are 
likely to come from the Planck satellite, we include spec- 
ifications for that satellite in Table [J We will also con- 
sider a cosmic variance limited experiment (Cosmic Var) 
to indicate the ultimate precision obtainable from CMB 
observations. In each case, we consider information from 
I < 2500, assuming that scales contaminated by sec- 
ondary anisotropics such as the Sunyaev-Zeldovich effect 
can be cleaned. We also exclude any information from 
B-mode polarization, although the B-modc signal gener- 
ated by weak gravitational lensing of E-modc polariza- 
tion does in principle contain information on neutrinos 
via the lensing potential power spectrum. 

Moving now to galaxy surveys, we may write the ap- 
propriate Fisher matrix as (45l . liq 



T^GAL 



ainP(k) 91nP(k)^^ ^ (P'k 

-Vcff(k); 



de. 



2(27r)= 



(5) 



where the derivatives are evaluated using the cosmologi- 
cal parameters of the fiducial model and V^s is the effec- 
tive volume of the survey, given by 



yoff(fc,/i) 



nP{k, fi) 



n(r)P(fc,M) 

-I 2 



1 



nP{k,fi) + 1 



(6) 



Here the galaxy survey is parametrized by the survey 
volume 14urvey and the galaxy density n(r), which in the 
last equality we assume to be uniform n within a red- 
shift bin. We assign an amplitude as,g for clustering of 
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galaxies on 8h Mpc"-'^ scales in each redshift bin. This is 
used to calculate the galaxy bias b = as^g/a^, where erg is 
the clustering of the dark matter. We take the observed 
power spectrum to be Pobs(fc) = b^P{k), neglecting pe- 
culiar velocity effects. Since the power spectrum is re- 
constructed from observations of galaxy positions using 
a fiducial cosmology P{k) aquires additional cosmologi- 
cal dependence arising from the Alcock-Paczynski effect. 
These help break some of the parameter degeneracies in 
the shape of the matter power spectrum by placing con- 
straints on H{z) and Dji{z). This is included in our 
analysis. 

In addition, we must specify fcmax, a cutoff on small 
scales to avoid the effects of non-linearity. We model this 
cut-off scale on the prescription of Seo and Eisenstein 
[4^, who assumed that non- linear effects become impor- 
tant on a scale fc„; = 7r/2/i?„; where i?„; is the scale 
on which averaged density fluctuations cr(i?„i) = 0.5. 
This represents a moderately conservative choice of cut- 
off. The analysis of Jeong and Komatsu [IJ suggests that 
using 3rd order perturbation theory the matter power 
spectrum may be accurately modeled well into the quasi- 
linear regime significantly extending the accessible range 
of fc-modes. More recently, Saito et al. [1^ have shown 
that, when the effect of neutrino masses on non- linear 
clustering is included in the calculation of the power spec- 
trum by 1-loop perturbation theory, the power spectrum 
suppression in the non-linear regime is enhanced. In this 
case, the combination of an enhanced signal and probing 
a larger volume by including the contribution of more 
scales increases the tightness of neutrino mass constraints 
by a factor of a few. Including these effects may allow 
for improvement over our results and indicates some of 
the complications of recovering cosmological information 
from small scales. Similar problems apply to the biasing 
of galaxies with respect to the background matter field, 
which may, on small scales, become scale dependent or 
be modified by the legacy of inhomogeneous reionization 
pollsoj. In this case, additional information from statis- 
tics beyond the power spectrum, e.g. the bispectrum 
[sH . [sa j. may be necessary to quantify and correct for 
these effects. 

Since much of the information about massive neutri- 
nos shows up as small scale modifications to the power 
spectrum our choice of the cutoff fcmax is important. Too 
small a value and we will throw away important informa- 
tion about massive neutrinos, too large a value and we 
will include scales affected by non-linear evolution and 
will underestimate our uncertainty. Note we will be con- 
sidering 21 cm information, which comes from smaller 
scales still, but since this comes from higher redshifts the 
non-linear scale is smaller. 

To indicate the spread of galaxy survey performance, 
we consider four galaxy surveys (with parameters mod- 
eled on those from [46[). Our basic galaxy survey, which 
we label SDSS, loosely corresponds to the SDSS LRG 
sample and is centered at z w 0.3. Next, we consider a 
larger galaxy survey centered at z = 1, labeled Gl. With 



TABLE II: Specification for galaxy surveys. For each galaxy 
survey, we specify the mean redshift 2, the survey volume 
Kiurvey, the mean galaxy density n in the redshift bin, the 
cutoff off scale fcmax, and the effective clustering of the galaxies 



Survey 


z 


Vsurvoy 


n 


^max 




SDSS 


0.3 


1.0 


10-^ 


0.1 


1.8 


Gl 


1.0 


1.8 


5 X lO-"* 


0.2 


1.0 


G2 


3.0 


0.5 


10'^ 


0.5 


1.0 


G3 


1.0 


40.0 


5 X 10"'' 


0.2 


1.0 



larger volume, higher galaxy density, and larger fcmax 
this survey represents a considerable improvement over 
SDSS. To illustrate the gain of going to higher redshift, 
we consider a Lyman break galaxy survey (G2) centered 
at z = 3. This has reduced volume, but compensates 
with higher galaxy density and larger fcmax- Finally, we 
consider an all-sky galaxy survey out to z = 1 (G3) ca- 
pable of surveying 40/i~'^Gpc^. This represents the limit 
of galaxy surveys in the local Universe and is enormously 
better than the other surveys that we consider. We list 
parameters for our galaxy surveys, along with our cho- 
sen cutoff scale, for each of our surveys in Table |TT1 We 
will consider surveys with spectroscopic redshifts, so that 
errors in galaxy redshifts are negligible. Photometric sur- 
veys may achieve the same level of constraint, although 
they require much larger volumes. The large redshifts er- 
rors greatly degrading sensitivity on small scales, which 
will hamper constraining neutrino masses. For example, 
the Large Scale Synoptic Telescope (LSST) is a full sky 
photometric survey of volume exceeding that of our G3, 
but since it relies upon photometric redshifts provides 
parameter constraints only a factor of ~ 3 better than 
Gl. 

Estimates of the error on a 21 cm power spectrum mea- 
surement are calculated using the Fisher matrix [H, [13] 

^ 47r2 aUk,^i) de, 80, ' 



In this equation, ^survey — D^AD{X'^/Ae) denotes the 
effective survey volume of our radio telescopes given the 
conformal distance D{z) to the center of the survey at 
redshift z, the depth of the survey AZ?, the observed 
wavelength A, and the effective collecting area of each 
antennae tile A^. 

The variance of a 21 cm power spectrum estimate for 
a single k-mode with line of sight component fc|| = /xfc, 
restricting ourselves to modes in the upper-half plane, in- 
cluding both sample variance and thermal detector noise, 
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and assuming Gaussian statistics, is given by |53l. l55l|: 

1 D'^AD ( A2 ~ 



1 



cld 



(8) 



The first term in this expression is the contribution from 
sample variance, while the second describes the thermal 
noise of the radio telescope. The thermal noise depends 
upon the system temperature Tgys, the survey bandwidth 
B, the total observing time <int- The effect of the config- 
uration of the antennae is encoded in the number density 
of baselines n_L (fc) that observe a mode with transverse 
wavenumber fcj, [53| . Observing a number of fields A'fioid 
further reduces the variance. 

In making predictions for 21 cm observations we will 
consider three separate instruments representing differ- 
ent levels of technological advancement. As an exper- 
iment representative of a 21 cm pathfinder experiment, 
we consider the Murtchison Widefield Array (MWA), cur- 
rently under construction in the Australian outback. For 
this we take 500 antennae distributed with a filled core 
of size 20m surrounded by the remainder of the antennae 
distributed with a r"^ profile out to 750m. To represent 
a mature 21 cm experiment, we consider an experiment 
modeled after the Square Kilometer Array (SKA). For 
this, we distribute 20% of a total of 5000 antennae within 
a 1km radius in the same way as the MWA -- a core 
surrounded by a distribution. These are then sur- 
rounded by a further 30% of the total antennae in a uni- 
form density annulus of outer radius 6km. The remain- 
der of the antennae are distributed at larger distances too 
sparsely to be useful for power spectrum measurements. 
Finally, we consider a hypothetical Fast Fourier Trans- 
form Telescope [56j . a future square kilometer collecting 
area array optimized for 21 cm observations. In contrast 
to MWA and SKA, which combine dipoles into tiles that 
are then cross-correlated, the FFTT performs an FFT 
to cross-correlate all the dipoles giving it considerably 
greater sensitivity. For the MWA and SKA, we assume 
that the effective collecting area of the instrument scales 
as A^, as expected for a dipole antennae. We summarize 
our choice of 21 cm experiment parameters in Table [Till 
For each experiment, we impose a cutoff at k — 2T:/y, 
where y is the conformal depth of the redshift slice, to 
account for the presence of foregrounds, whose removal 
will reduce the information available on large scales [5^ . 
We also impose fcmax = 3hMpc~^ as the non- linear scale 
and exclude larger k values from the analysis. 

To give a detailed comparison of these experiments, 
and so better understand why they result in different 
parameter constraints, it is useful to define and plot the 
effective volume observed. It has previously been pointed 
out that 21 cm experiments potentially explore a larger 
volume [57} . The effective volume is defined explicitly for 
galaxy surveys in Equation ^ and we define an equiva- 
lent quantity for 21 cm experiments (essentially weight- 
ing the survey volume by the signal to noise) by analogy 



TABLE III: Low-frequency radio telescopes and their param- 
eters. For MWA we assume a single redshift slice centered at 
z = 8, while for SKA and FFTT we divide the full redshift 
range into five redshift slices of thickness Az ~ 0.5. 



Array 




Atot 


-^min 




B 


Tint 


z 






(lO^m^) 


(m) 


(km) 


(MHz) 


(hr) 




MWA 


500 


7.0 


4 


1.5 


8 


4000 


7.8-8.2 


SKA 


5000 


600 


10 


5 


8 


4000 


7.8-10.3 


FFTT 


10^ 


10^ 


1 


1 


8 


4000 


7.8-10.3 



to Equation ^ 



21cm 
off 



{k,n) 



Pi 



Tl{k,n) 



14 



survey • 



(9) 



The different effective areas for the experiments that we 
consider are plotted in Figure [31 Looking first at the 
galaxy surveys, we see that Gl probes more volume on 
all scales than SDSS leading to improved parameter con- 
straints. G2 is less sensitive on larger scales, but is able 
to probe deeper than SDSS leading to extra constraining 
power on parameters, like M^, which affect small scales. 
MWA's concentrated antennae distribution gives it good 
sensitivity to modes fc ~ 0.1 Mpc~^, but drops in sen- 
sitivity dramatically on small scales, where it is noise 
dominated. SKA's greater collecting area and more dis- 
persed layout allows it to constrain scales fc 1 Mpc~^. 
The FFTT has a much larger effective volume than any 
of these experiments and is effectively sample variance 
limited out to fc 1 Mpc"""^. This tremendous volume 
gives it the power to constrain cosmological parameters 
at unprecedented levels. It is instructive to compare it 
with G3, which probes similarly large volumes, but whose 
sensitivity drops rapidly on scales fc > 0.1 Mpc~^. This 
different behavior makes all the difference when it comes 
to constraining M^. One unfortunate feature is that fore- 
ground removal, which removes the large scales, is likely 
to prevent a 21 cm experiment observing the imprint of 
neutrinos on the power spectrum at both large scales 
(fc ^ 0.001ft. Mpc~^), where there is little effect, and 
small scales (fc > OAh Mpc~^), where the power spec- 
trum is highly suppressed. 

Why is there this difference in sensitivity on small 
scales? For a galaxy survey to have useful sensitivity it 
is important that there be a large enough density of tar- 
get galaxies, i.e. we require P{k)n > 1. If the observed 
galaxies are sparsely distributed within the survey vol- 
ume, then they only poorly trace small scale features in 
the dark matter power spectrum. Consequently, galaxy 
surveys are at the mercy of the details of galaxy forma- 
tion and lose sensitivity once pushed to scales far below 
the typical galaxy spacing. In contrast, the sensitivity of 
21 cm experiments on small scales is a strong function of 
the distribution of antennae - a design feature under hu- 
man control. If the array lacks a sufficiently high density 
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FIG. 3: Comparison of angle averaged effective volume for 
galaxy and 21 cm surveys. SDSS (thin solid curve), Gl (thin 
dotted curve), G2 (thin short dashed curve), G3 (thin long 
dashed curve), MWA (thick solid curve), SKA (thick dotted 
curve), and FFTT (thick dashed curve). 



of long baselines it becomes noise dominated at large k 
values. By close filling the array to large baselines, as in 
FFTT, sensitivity of a 21 cm experiment can be large for 
even very large k values. Since hydrogen gas is present 
on all scales there is no need to worry about an equiva- 
lent shot noise component. While the small scale 21 cm 
power spectrum might be modified on small scales by as- 
trophysical effects, given the high sensitivity achievable 
it is worth exploring the use of 21 cm studies further. 



V. RESULTS 
A. Total neutrino mass 

Table ITVl summarizes the Fisher matrix parameter con- 
straints for galaxy surveys and 21 cm experiments alone 
and in combination with Planck. We also list results for 
a CosmicVar and FFTT. Our parameter constraints on 
Ml, are consistent with the predictions of ,37j for Planck 
and SDSS. 

As we have shown 21 cm experiments benefit from ac- 
cessing a larger effective volume than most galaxy sur- 
veys. This is mostly a consequence of the well known 
fact that there is more volume available at higher red- 
shift. The MWA performs comparatively to large galaxy 
surveys such as Gl and G2. With the SKA we begin to 
exploit the benefit of larger volume seeing the neutrino 



FIG. 4: 2a constraints on M,^ (solid curve) from 
SKA+PLANCK (dotted curve) and FFTT+PLANCK 
(dashed curve) for the normal (thin curves) and inverted 
(thick curves) hierarchies. 



constraint halve to ctm^ = 0.16. At this level, SKA, by 
itself, could detect our fidicial neutrino mass M^, = 0.3, 
although barely at the 2 — cr level. Adding in Planck to 
MWA and SKA leads to substantial improvement and 
allowing a, 4 — a detection of neutrino mass. The FFTT's 
high sensitivity allows it to constrain neutrino mass at 
the level of (Tm^ = 0.01 by itself. This is sufficient to de- 
tect either normal or inverted mass hierarchies and would 
open up an era of precision cosmological measurements of 
neutrino mass. In general, combining Planck with other 
21 cm experiments helps with M^, constraints, although 
this is a very weak effect for FFTT. Planck is also vital 
for breaking the degeneracy between bias and As- 

Figure [3] illustrates the effect of changing on the 
parameter constraints. Since only SKA and FFTT are 
capable of placing useful constraints on M^, we restrict 
to plotting curves for these experiments. In general, the 
mass constraints are relaxed as we move to smaller M^, 
where the impact of neutrino mass is lessened, but the 
constraint does not depend upon the hierarchy. 

Based upon these results, SKA-|-Planck should be ca- 
pable of detecting neutrino mass at the 2 — cr mass pro- 
vided that > 0.15 eV. FFTT+Planck is capable of 
detecting neutrino mass for all masses above the mini- 
mum mass in both normal and inverted hierarchies. Fur- 
ther, FFTT-|-Planck has sufficient sensitivity to be sensi- 
tive to the hierarchy itself. From Figure [H we see that it 
could measure masses in the range M^, < 0.11 eV, below 
the minimum mass allowed for the inverted hierarchy and 
so confirm the existence of a normal hierarchy with 
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TABLE IV: Fiducial parameter values and 1 — a experimental uncertainties for basic 8 parameter model. Dashes indicate 
parameters not relevant for that experiment; oo indicates parameters that are relevant, but not constrained. 















r 






bsDSS 


Tt{8) Tb{8.5) Tb{9) Tt{9.5) Tt{10) 


Fiducial 


0.15 


0.023 


0.7 


-1 0.95 


27 


0.1 


0.24 


0.3 


2.9 


26 


26 


27 


27 


28 


SDSS 


0.154 


0.0363 


0.0548 


- 0.425 


oo 


- 


- 


2.54 


oo 












Gl 


0.0294 


0.00722 


0.00988 


- 0.07 


oo 


- 


- 


0.431 


oo 












G2 


0.0216 


0.0058 


0.00648 


- 0.0321 


oo 


- 


- 


0.371 


oo 












G3 


0.00623 


0.00153 


0.0021 


- 0.015 


oo 


- 


- 


0.091 


oo 












MWA 


0.0279 


0.0075 


0.0139 


- 0.046 


oo 


- 


- 


0.57 




oo 










SKA 


0.00043 


0.00056 


0.0019 


- 0.005 


oo 


- 


- 


0.16 




oo 


oo 


oo 


oo 


oo 


FFTT 


3.8e-05 


4.77e-05 


9.18e-05 


- 0.0002 


oo 


- 


- 


0.0089 




oo 


oo 


oo 


oo 


oo 


Planck 


0.0045 


0.00022 


0.04 


- 0.0071 


0.26 


0.0048 


0.01 


0.38 




+SDSS 


0.00236 


0.000204 


0.0216 


- 0.0068 


0.25 


0.00448 


0.01 


0.217 


0.197 












+G1 


0.000924 


0.000196 


0.00749 


- 0.0066 


0.243 


0.00439 


0.00989 


0.101 


0.067 












+G2 


0.00081 


0.0002 


0.0063 


- 0.0063 


0.24 


0.0044 


0.0092 


0.10 


0.12 












+G3 


0.00037 


0.00016 


0.0016 


- 0.0051 


0.24 


0.0043 


0.0081 


0.046 


0.029 












+MWA 


0.0013 


0.0002 


0.011 


- 0.0063 


0.24 


0.0044 


0.0095 


0.13 




0.74 










+SKA 


0.00029 


0.00014 


0.0016 


- 0.003 


0.23 


0.0043 


0.0041 


0.075 




0.46 


0.47 


0.48 


0.49 


0.49 


+FFTT 


3.72e-05 


3.58e-05 


9.13e-05 


- 0.0002 


0.23 


0.0043 


0.003 


0.0075 




0.12 


0.12 


0.12 


0.13 


0.13 


Cosmic Var 


0.00148 


4.12e-05 


0.014 


- 0.0024 


0.12 


0.0021 


0.0027 


0.14 














+FFTT 


3.28e-05 


2.26e-05 


8.86e-05 


- 0.0002 


0.11 


0.0021 


0.00096 


0.0059 




0.06 


0.06 


0.07 


0.07 


0.07 



in this range. 

Table fVl summarizes the Fisher matrix parameter con- 
straints for galaxy surveys and 21 cm experiments alone 
and in combination with Planck when the equation of 
state w is allowed to vary. This significantly weakens the 
constraints from galaxy surveys, but has little effect on 
the neutrino constraints for 21 cm experiments. Since 
21 cm experiments are generally well within the mat- 
ter dominated regime they are only sensitive to the dark 
energy through the angular diameter distance. This is 
essentially the same measurement as made by the CMB. 
Once Planck is added to galaxy or 21 cm surveys very 
similar constraints as in the case of ACDM are ob- 
tained. This extra data set is very powerful at breaking 
the degeneracy. 

We summarize these result in Figure [5l where we plot 
68 — % confidence intervals in the w — M^ plane. The con- 
tours for galaxy surveys show a clear diagonal alignment 
indicating significant degeneracy between w and Mi,. 
Hence better measurements of either quantity translate 
into better measurements of the other. In contrast, the 
contours for 21 cm experiments are less closely aligned 
with the axes indicating little degeneracy. This is as ex- 
pected since the 21 cm experiments are focused at high 
redshifts where dark energy has little effect. This also 
explains why 21 cm experiments are generally quite poor 
at constraining w. Note that since these observations 
take place at very different redshifts the axes of the er- 
ror ellipse are not aligned implying that combining all 
data sets would further reduce the maximum constraint 
achievable. 



Finally, it is worth noting that constraints on Ifjc may 
themselves be interesting from the point of view of con- 
straining neutrino physics. In this analysis, as in most, 
we have treated Yho as a free parameter to be varied in- 
dependently. In the framework of cosmology, however, 
once fifc/i^ is known predictions for Yho can be made us- 
ing the machinery of Big Bang Nucleosynthesis (BBN). 
A discrepancy between the prediction and observation 
would point to an early modification of general relativity 
or, more interesting for neutrinos, the presence of ex- 
tra relativistic degrees of freedom iV^.e// during BBN, 
such as a sterile neutrino. The power of this additional 
constraint was pointed out by [58| and here we simply 
comment on the applicability of 21 cm observations. Cur- 
rent constraints on Ihc come from observations of Helium 
lines in metal poor HII regions. By extrapolating the ob- 
served abundances back to zero metallicity a primordial 
abundance can be estimated. These observations place 
bounds Yue = 0.2477 ± 0.0029 [59], but there is consid- 
erable systematic uncertainty relating to, amongst other 
problems, the unknown temperature profile of the HII 
region. Given these systematic effects it is unclear that 
these constraints can be improved much beyond an addi- 
tional factor of two further. Additional constraints from 
more reliable sources are clearly important. Table IIVI 
shows that Planck only weakly constrains Yhb , since Yne 
modifies the power spectrum primarily on small scales 
where Silk damping reduces the amplitude of the signal. 
Adding additional information on small scales through 
21 cm observations helps break this degeneracy leading 
to factor of 2 (3) improvements for SKA (FFTT). The 
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TABLE V: Fiducial parameter values and 1 — a experimental uncertainties for basic 8 parameter model+variation in w. Dashes 
indicate parameters not relevant for that experiment; oo indicates parameters that are relevant, but not constrained. 











w 


Us 




r 




M, 


b.SDSS 


Tt{8) Tt{8.5) Tt{9) Tt{9.5) Tt{10) 


Fiducial 


0.147 


0.023 


0.7 


-1 


0.95 


26.6 


0.1 


0.24 


0.3 


2.3 


26 


26 


27 


27 


28 


SDSS 


0.456 


0.083 


0.117 


1.21 


0.503 


CO 


- 


- 


6.16 


oo 












Gl 


0.119 


0.0207 


0.0358 


0.574 


0.174 


oo 


- 


- 


1.28 


oo 












G2 


0.0354 


0.00593 


0.295 


1.22 


0.0482 


oo 


- 


- 


1.01 


oo 












G3 


0.0252 


0.00438 


0.0076 


0.122 


0.037 


oo 


- 


- 


0.272 


oo 












MWA 


0.0317 


0.00761 


0.972 


3.13 


0.0487 


oo 


- 


- 


0.749 




oo 










SKA 


0.00191 


0.00056 


0.234 


0.747 


0.0054 


oo 


- 


- 


0.175 




oo 


oo 


oo 


oo 


oo 


FFTT 


4.94e-05 


4.77e-05 


0.0045 


0.014 


0.0002 


oo 


- 


- 


0.009 




oo 


oo 


oo 


oo 


oo 


Planck 


0.0045 


0.00024 


0.068 


0.18 


0.0074 


0.26 


0.0048 


0.011 


0.38 




+SDSS 


0.0033 


0.00024 


0.023 


0.11 


0.0074 


0.254 


0.0046 


0.0103 


0.272 


0.197 












+G1 


0.0016 


0.00021 


0.013 


0.081 


0.0068 


0.245 


0.0044 


0.01 


0.136 


0.067 












+G2 


0.00089 


0.00022 


0.037 


0.149 


0.0067 


0.243 


0.0044 


0.0099 


0.104 


0.16 












+G3 


0.00051 


0.00016 


0.003 


0.021 


0.0051 


0.24 


0.0043 


0.0081 


0.052 


0.029 












+MWA 


0.00146 


0.00021 


0.053 


0.17 


0.0066 


0.242 


0.0044 


0.0101 


0.144 




1.45 










+SKA 


0.00029 


0.00014 


0.020 


0.065 


0.003 


0.236 


0.0043 


0.0044 


0.080 




0.71 


0.72 


0.74 


0.75 


0.76 


+FFTT 


4.23e-05 


3.97e-05 


0.004 


0.011 


0.0002 


0.23 


0.0043 


0.0030 


0.0075 




0.14 


0.15 


0.15 


0.15 


0.16 


Cosmic Var 


0.00244 


4.16e-05 


0.030 


0.033 


0.0024 


0.124 


0.0023 


0.0028 


0.222 














+FFTT 


3.3e-05 


2.58e-05 


0.0024 


0.0076 


0.0002 


0.111 


0.0021 


0.0011 


0.0068 




0.09 


0.09 


0.1 


0.1 


0.1 



constraint from Planck + FFTT of 0.0028 is comparable 
with current observational bounds. Note that although, 
through the baryon component, Yhc affects the matter 
power spectrum on small scales the effect is very small. 
We find that by themselves SKA (FFTT) can place con- 
straints on Yhc of 0.17 (0.014). This extra information 
offers only a small improvement over that contained in 
the CMB. 



B. Neutrino mass splitting 

We have seen in the previous subsection that SKA, G3, 
and FFTT all possess the sensitivity needed to provide 
a detection of neutrino mass at high sensitivity (at least 
for some range of Mi,), ft is interesting to ask if any 
of these experiments can go further and begin to con- 
strain the masses of individual neutrino mass states. As 
mentioned in i jllll the power spectrum is primarily sensi- 
tive to the total neutrino energy density, but is modified 
at the 0.2% level by mass differences between the three 
neutrino states. Precise power spectrum measurements 
could, in principle, make this measurement. This would 
conclusively identify neutrinos as acting in a cosmological 
context. 

We begin our analysis by plotting the 1 — cr power spec- 
trum errors for G3, SKA, and FFTT, which are shown 
in Figure [6] In each case, we consider the errors on a 
fiducial model with M^, = 0.3 distributed in the normal 
hierarchy (dashed line). For comparison, we plot curves 
for the power spectrum in the case of the inverted hierar- 



chy (dotted curve) and for three neutrinos with degener- 
ate masses (solid curve). These curves give an indication 
of how precise power spectrum measurements must be to 
have a hope of distinguishing individual neutrino masses. 
In the top panel of Figure [SI we plot errors for G3 on 
the power spectrum at z = 1. It is clear that while G3 
measures the power spectrum at the 0.1% level it is still 
not precise enough to distinguish the different mass split- 
tings. This is consistent with the findings of Lesgourgues 
et al. [IJ- Note that non-linearities in galaxy bias and 
the underlying density field would interfere with measure- 
ments of the power spectrum on scales k > 0.2h Mpc~^ 
anyway, further restricting the ability of a low redshift 
galaxy survey to constrain neutrino mass splittings. 

In the bottom panel of Figure [51 we consider power 
spectrum errors from SKA and FFTT at z = 8. First 
note that the differences between the three different mod- 
els has different scale dependence than at z = 1. Unfor- 
tunately, this does not greatly increase the chances of 
detecting the mass splitting. Next, since the non-linear 
scale at z = 8 is /c > 3ft. Mpc~^ the whole of the plot- 
ted region is potentially accessible to experiment with- 
out having to worry about complications of a non-linear 
density field. We will discuss complications arising from 
small scale modifications due to inhomogeneous radia- 
tion fields to the 21 cm power spectrum later. SKA 
has almost the sensitivity needed to differentiate between 
different mass splittings at the 1 — a level on scales 
0.4 < A: < Ih Mpc-\ The FFTT does visibly better, 
having error bars smaller than the difference in the three 
power spectra on scales k > 0.2ft, Mpc~^. This moti- 
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FIG. 5: 68-% confidence ellipses in the — w plane from 
PLANCK in combination with SDSS (thin dotted), G2 (thin 
short dashed), Gl (thin solid), G3 (thin long dashed), MWA 
(thick dotted), SKA (thick dashed), and FFTT (thick solid). 



vates the next part of our analysis, where we analyze the 
Fisher matrix constraints on individual neutrino masses 
for these experiments, since degeneracies between differ- 
ent parameters affect the detectability of neutrino mass. 



FIG. 6: Top panel: 1-a fractional errors on the power spec- 
trum at z = 1 for G3. Bottom panel: 1-a fractional errors 
on the power spectrum at 2 = 8 for SKA (thin black points), 
and FFTT (thick black points). In each panel, we plot the 
ratio of P{k) for three degenerate masses (dotted curve) or 
inverted hierarchy (solid curve) to P{k) for the normal hi- 
erarchy neutrinos (as in Figure [2} and = 0.3 eV. Only 
the FFTT has the sensitivity to detect the mass splittings at 
a significant level. Nonlinearities are expected to affect the 
z = 1 power spectrum on scales smaller than k > 0.2h Mpc~^ 
and the z = 8 power spectrum on scales k > 3h Mpc~^. 



TABLE VI: Neutrino mass 1 — a uncertainties. These uncer- 
tainties are calculated using the full 7 parameter plus three in- 
dependent neutrino masses Fisher matrix and then projected 
onto the reduced parameter set listed. 





mi m2 ma 


mi -1- m2 ms 




Fiducial 


0.081 0.09 0.13 


0.17 0.13 


0.3 


Planck 
+G3 
-fSKA 
+FFTT 


0.922 0.950 0.552 
0.654 0.78 0.295 
0.645 0.631 0.405 
0.037 0.03 0.03 


0.633 0.550 
0.296 0.265 
0.41 0.405 
0.03 0.029 


0.44 
0.046 
0.077 
0.0076 


Cosmic Variance 
+G3 
+SKA 
_I_FFTT 


0.285 0.282 0.209 
0.247 0.272 0.139 
0.265 0.264 0.194 
0.036 0.03 0.03 


0.235 0.209 
0.15 0.137 
0.192 0.194 
0.029 0.028 


0.157 
0.027 
0.037 
0.006 



We now examine in detail the ability of our fiducial ex- 
periments to constrain individual neutrino masses. We 
take as our fiducial model three neutrino mass states dis- 
tributed according to the normal hierarchy with M^^ — 
0.3 eV. This fiducial mass is sufficiently large that a mea- 
surement of Ml, alone would not allow either hierarchy 



to be ruled out. Only if individual neutrino masses could 
be measured could the hierarchy be established. 

While CAMB has sufficient precision to calculate the 
effect on the power spectrum of individual neutrino 
masses, the resulting differences in the power spectrum 
are only slightly larger than numerical uncertainties in 
the code (just visible as numerical noise in Figure [2]). 
We address this by running CAMB at its highest accu- 
racy settings and visually inspecting the power spectrum 
derivatives to ensure that there are no distinctive glitches. 
Large numerical noise could easily break the degeneracy 
between different neutrino masses yielding erroneously 
small error bars. As a numerical check, by projecting 
from errors on (mi, TO2, 7713) onto M^, we see that these 
results are consistent with the results of the previous sec- 
tion. 

Parameter constraints are shown in Table IVII for indi- 
vidual neutrino masses. From the fisher matrix for three 
individual neutrino masses, we project onto the {mi+m2, 
7713} and {M,y} parameter spaces. Since the mass split- 
ting between mi and m2 is < 0.01 eV there exists a 
strong degeneracy between these two masses for most 
experiments. In contrast the sum of these two nearly 
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degenerate masses is much better constrained. By con- 
sidering the sum, we get a better handle on whether the 
experiments could distinguish between normal and in- 
verted hierarchy by isolating whether one or two neutrino 
mass states lie above the largest mass splitting. 

Neither G3 or SKA is able to place significant con- 
straints on the individual neutrino masses. FFTT is 
able to distinguish from zero all three mass states at the 
> 3 — (7 level. We conclude that an FFTT type 21 cm ex- 
periment in conjunction with at least Planck level CMB 
data might someday make a concrete detection of indi- 
vidual neutrino masses in the cosmological context. 

For comparison with these results, the inverted 
hierarchy with AI^ = 0.3 eV satisfies (mi, m2, 
ra3)=(0.066,0.113,0.122) eV differing from the normal hi- 
erarchy by (Ami, Ato2, Ato3)=(0.015,0.023,0.008) eV. 
These differences are too small for any of the experiments 
to distinguish between the two hierarchies (see Figure 
[7]). Thus although individual neutrino masses may be 
measured, the precision required to distinguish normal 
and inverted hierarchies directly at the level of individ- 
ual masses is difficult to achieve. It would be possible in 
principle by increasing the collecting area or integration 
time of the FFTT, although the design we consider is 
already quite futuristic. 



VI. CONCLUSIONS 

For the foreseeable future cosmology is likely to provide 
insights into neutrino masses, constraining the underly- 
ing particle physics. In this paper, we have explored the 
possibility of using low-frequency radio observations of 
the 21 cm line of hydrogen during the epoch of reioniza- 
tion to constrain neutrino masses and compared with fu- 
ture galaxy surveys. In general, we find that the outlook 
for 21 cm experiments is promising since these experi- 
ments are capable of probing small scale modifications of 
the matter power spectrum with higher sensitivity than 
galaxy surveys. This is in part a consequence of going 
to higher redshifts where non-linear modifications to the 
matter power spectrum are expected to be smaller. It is 
also a consequence of the sensitivity of 21 cm experiments 
being dominated on small scales by the distribution of 
antennae. If the radio array is constructed with a large 
enough density of long base lines it can be made sensi- 
tive to large fc-modes. In contrast, galaxy surveys are 
reliant upon there being a large enough density of target 
galaxies so that power spectrum measurements are not 
dominated by shot noise at large k values. 

We find that 21 cm experiments such as SKA should 
be competitive with galaxy surveys for constraining the 
total neutrino mass. In combination with Planck, SKA 
should be able to measure the effects of neutrino mass 
at the 2 — (T level if > 0.15. These constraints are 
largely independent of uncertainties in the dark energy, 
which affects lower redshift probes of the matter power 
spectrum. In the further future, an instrument with a 




FIG. 7: 68% confidence contours on individual neutrino mass 
states nil, m2, and ms. We show contours for Cosmic Var in 
combination with G3 (long dashed curve), SKA (short dashed 
curve), and FFTT (solid curve). Fiducial neutrino masses 
are indicated for the normal hierarchy (diagonal crosses), 
also shown are the values for the inverted hierarchy (verti- 
cal crosses) for = 0.3 eV. Cosmic Var-|-FFTT is able to 
make a weak detection of individual neutrino masses. Insert 
in the bottom right hand corner shows the error ellipse in 
the mi + m2-m2. parameter space and illustrates the strong 
degeneracy hampers placing tight constraints on individual 
masses. 



square kilometer of collecting area able to cross-correlate 
all its antennae, such as FFTT, could measure Mi, at the 
better than 2 — a for all values of M^, currently allowed by 
the observed mass splittings. This is exciting as it raises 
the possibility of using cosmology to distinguish between 
the normal and inverted neutrino mass hierarchies. Pre- 
vious work [ST^I had shown that even an all-sky galaxy 
survey out to 2 = 1 is unlikely to be able to measure 
individual neutrino masses. We obtain a similar conclu- 
sion for a 21 cm experiment such as SKA, but note that 
a successor instrument, such as FFTT, might indeed be 
able to distinguish individual masses. Our calculations 
show that FFTT could make > 2 — cr detections of each 
of three neutrino mass states. This would be valuable 
for making a measurement of neutrino masses, if labora- 
tory experiments have not already done so, or confirming 
that we are indeed seeing the effect of neutrinos in the 
cosmological context. 

Whether 21 cm experiments live up to the promise 
shown by this work depends on a number of factors. First 
of all is the very real problem of astrophysical foregrounds 
[boI . which must be overcome before 21 cm fluctua- 
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tions can be measured and used for cosmolog;^ The out- 
look for foreground removal is positive [63 , iGSj and will 
be better understood once the first generation of instru- 
ments, e.g. LOFAR, MWA, 21CMA, are completed and 
collect data. Next is the possibility that modifications 
to the 21 cm power spectrum from rcionization or from 
astrophysical radiation fields degrade the accuracy with 
which the underlying density field can be measured. We 
have taken the optimistic point of view that this will not 
affect parameter constraints in order to get a sense of 
how well one might hope to do. In light of recent work 
by Mao et al. [32] , it is likely that these effects will be an 
issue, but not a show stopper. Mao et al. [s^ find that 
needing to fit for unknown parameters describing rcion- 
ization degrades cosmological parameter constraints by 
between 10 — 50%. This would not prevent the detec- 
tion of individual neutrino masses that we have described 
here, although clearly it would reduce their significance. 
In summary, it seems possible that future 21 cm experi- 
ments will provide highly complementary constraints on 
neutrino masses and could push beyond those obtained 
from galaxy surveys. There are clearly challenges ahead 



for this technique, but if overcome, the future of neutrino 
cosmology may be bright. 
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